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a: 

<-j ! The 0(a) electroweak radiative corrections to the process — > — > £ ± f 

= e, /x) are calculated. The C(a) corrections can be decomposed into 
separately gauge invariant contributions to the W boson production and decay 
$_i 1 processes. Factorizing the collinear singularity associated with initial state 



photon radiation into the parton distribution functions, we find that initial 
state corrections have a significantly smaller effect than final state radiative 
corrections. We study in detail the effect of electroweak radiative corrections 
on a number of interesting observables: the W transverse mass distribution, 
the W to Z transverse mass ratio, the charge asymmetry of leptons in W — > Iv 
decays, as well as the W production cross section and the W to Z cross section 
ratio. We also investigate how experimental lepton identification requirements 
change the effect of the electroweak corrections. 
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I. INTRODUCTION 



The Standard Model of electroweak interactions (SM) so far has met all experimental 
challenges and is now tested at the 0.1% level |IJ. However, there is little direct experimental 
information on the mechanism which generates the masses of the weak gauge bosons. In the 
SM, spontaneous symmetry breaking is responsible for mass generation. The existence of 
a Higgs boson is a direct consequence of this mechanism. At present the negative result of 
direct searches performed at LEP2 imposes a lower bound of Mh > 87.6 GeV on the Higgs 
boson mass. Indirect information on the mass of the Higgs boson can be extracted from 
the M H dependence of radiative corrections to the W boson mass, M w , and the effective 
weak mixing angle, sin 2 O^f*. Assuming the SM to be valid, a global x 2_ fit to all available 
electroweak precision data yields a (one-sided) 95% confidence level (CL) upper limit on Mh 
of 408 GeV @. 

Future more precise measurements of Mw and the top quark mass, m top , will lead to more 
accurate information on the Higgs boson mass 0-0. Currently, the W boson mass is known 
to ±65 MeV from direct measurements, whereas the uncertainty of the top quark mass is 
±5.2 GeV §. With a precision of 30 MeV (10 MeV) for the W mass, and 2 GeV for the top 
quark mass, Mh can be predicted from a global analysis with an uncertainty of about 30% 
(15%) Comparison of these indirect constraints on Mh with the results from direct 

Higgs boson searches at LEP2, the Tevatron collider, and the Large Hadron Collider (LHC) 
will be an important test of the SM. They will also provide restrictions on the parameters 
of the Minimal Supersymmetric extension of the Standard Model (MSSM) ||. 

A significant improvement in the W mass uncertainty is expected in the near future from 
measurements at LEP2 [10] and the Fermilab Tevatron pp collider ||. The ultimate precision 
expected for My/ from the combined LEP2 experiments is approximately 40 MeV ||10|| . At 
the Tevatron, integrated luminosities of order 1 fb _1 are envisioned in the Main Injector 
Era (Run II), and one expects to measure the W mass with a precision of approximately 
50 MeV H per experiment. The prospects for a precise measurement of Mw would further 
improve if a significant upgrade in luminosity beyond the goal of the Main Injector could be 



realized. With recent advances in accelerator technology ITT], Tevatron collider luminosities 
of order 10 33 cm~ 2 s _1 may become a reality, resulting in integrated luminosities of up to 
10 fb _1 per year. With a total integrated luminosity of 30 fb _1 , one can target a precision 
of the W mass of 15 - 20 MeV . A similar or better accuracy may also be reached at the 
LHC 0]. 

In order to measure the W boson mass with high precision in a hadron collider envi- 
ronment, it is necessary to fully understand and control higher order QCD and electroweak 
(EW) corrections. A complete calculation of the full 0(a) electroweak radiative corrections 
to pp —>■ W ± — > i^v (£ = e, /i) has not been carried out yet. In a previous calculation, 
only the final state photonic corrections were included |I1||I4]| , using an approximation in 
which the sum of the soft and virtual parts is indirectly estimated from the inclusive 0(a 2 ) 
W —>■ iv{^{) width and the hard photon bremsstrahlung contribution. The unknown part of 
the 0(a) electroweak radiative corrections, combined with effects of multiple photon emis- 
sion (higher order corrections), have been estimated to contribute a systematic uncertainty 



of SM\v = 15 — 20 MeV to the measurement of the W mass 15,16 



In this paper we present a new and more accurate calculation of the 0{a) EW correc- 
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tions to resonant W boson production in hadronic collisions. Our calculation is based on 
the full set of 0(a 3 ) Feynman diagrams, and includes both initial and final state radiative 
corrections, as well as the contributions from their interference. Final state charged lepton 
mass effects are included in the following approximation. The lepton mass regularizes the 
collinear singularity associated with final state photon radiation. The associated mass sin- 
gular logarithms of the form In (s/mf), where s is the squared parton center of mass energy 
and me is the charged lepton mass, are included in our calculation, but the very small terms 
of 0{m\/s) are neglected. 

To perform our calculation, we use a Monte Carlo method for next-to-leading-order 
(NLO) calculations similar to that described in Ref. With the Monte Carlo method, it 
is easy to calculate a variety of observables simultaneously and to simulate detector response. 
Calculating the EW radiative corrections to resonant W boson production, the problem 
arises how an unstable charged gauge boson can be treated consistently in the framework of 
perturbation theory. This problem has been studied in Ref. JT^] with particular emphasis 
on finding a gauge invariant decomposition of the EW 0(a) corrections into a QED-like 
and a modified weak part. Unlike the Z boson case, the Feynman diagrams which involve a 
virtual photon do not represent a gauge invariant subset. In Ref. fll8| , it was demonstrated 
how gauge invariant contributions that contain the infrared (IR) singular terms can be 
extracted from the virtual photonic corrections. These contributions can be combined with 
the also IR-singular real photon corrections in the soft photon region to form IR-finite gauge 
invariant QED-like contributions corresponding to initial state, final state and interference 
corrections. The collinear singularities associated with initial state photon radiation can 
be removed by universal collinear counter terms generated by "renormalizing" the parton 
distribution functions |19|J20|| , in complete analogy to gluon emission in QCD. A similar 
strategy has been employed in a recent calculation of the 0(a) QED corrections to Z boson 



production in hadronic collisions [21 



The technical details of our calculation are described in Sec. II. We first extract the 
collinear behaviour of the partonic cross section for both the initial and the final state 
corrections. Then, we define the quark distribution functions in next-to-leading order QED 
within the QED DIS and the QED MS factorization scheme when a finite quark mass 
is used to regulate the collinear singularities. Finally, we provide explicit formulae for the 
0(a 3 ) differential cross section and perform various consistency checks. 

Numerical results for pp collisions at ^/s = 1.8 TeV are presented in Sec. III. In hadron 
collider experiments, W bosons are identified by their leptonic decays, W — > tv. Since the 
neutrino escapes undetected, the Iv invariant mass cannot be reconstructed, and one must 
resort to other kinematic variables for the measurement of Myy. The observable which cur- 
rently provides the best measurement of Mw is the distribution of the transverse mass, Mp. 



The Mt distribution sharply peaks at Mw, and is rather insensitive to QCD corrections |22 



Alternative measurements of the W mass are provided [16j by the lepton p? distribution 
which peaks at M w /2, and the W/Z transverse mass ratio |23| , |2^| . Due to the mass singular 
logarithms associated with final state photon bremsstrahlung in the limit where the photon 
is emitted collinear with the charged lepton, the distributions which are sensitive to M w , the 
W production cross section and the W to Z cross section ratio in presence of cuts, and the 
charge asymmetry of leptons in W decays are significantly affected by the 0(a) electroweak 
radiative corrections. 
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The size of the radiative corrections strongly depends on the detector resolution. In 
Sec. Ill, using a simplified model of the D0 detector as an example, we also investigate how 
the finite resolution of realistic detectors affects the electroweak radiative corrections. Elec- 
trons and photons which are almost collinear are difficult to discriminate, and the momenta 
of the two particles are thus recombined into an effective electron momentum [f[|,|16|] if their 
separation in the pseudo-rapidity - azimuthal angle plane is below a critical value. This 
procedure completely eliminates the mass singular logarithms and, therefore, strongly re- 
duces the size of the 0(a) corrections. In contrast, photons which are almost collinear with 
muons are rejected if they are too energetic P"o] , |2o| which results in residual mass singular 
logarithmic corrections to observable quantities in W production. Finally, our conclusions 
are presented in Sec. IV. 



II. O(o) ELECTROWEAK CORRECTIONS TO W PRODUCTION IN 

HADRONIC COLLISIONS 

The calculation presented here employs a combination of analytic and Monte Carlo 
integration techniques^ Details of the method can be found in Ref. [fT7|| . The Feynman 



diagrams contributing to W boson production in hadronic collisions to 0(a 3 ), 

Qi(pi)qAPi') - ^ + (?)(7) - MPf)e + (Pf')(l(k)) 

are shown in Fig. 0. Since we are interested in the cross sections in the vicinity of the W 
resonance, the W, Z box diagrams can be neglected as non-resonant contributions of higher 
order in perturbation theory, and thus are not depicted in Fig. |l|. The calculation of Iv 
production in hadronic collisions at 0(a 3 ) includes contributions from the square of the 
Born graphs, the interference between the Born diagrams and the virtual one loop graphs, 
and the square of the real emission diagrams. 

Our treatment of the 0(a) corrections to W boson production in the resonance region is 
based on the calculation presented in Ref. fl8|], which we outline here. Unlike the Z boson 



case, the Feynman diagrams of Fig. [I] which involve a virtual photon do not represent a gauge 



invariant subset. In Ref. iTSfl , it was demonstrated how gauge invariant contributions that 
contain the infrared singular terms can be extracted from the virtual photonic corrections. 
These contributions can be combined with the also IR-singular real photon corrections in the 
soft photon region to form IR-finite gauge invariant QED-like contributions corresponding 
to initial state, final state and interference corrections. The soft photon region is defined by 
requiring that the photon energy in the parton center of mass frame, E 7 , is _E 7 < E cut = 
5 s y~§/2. In this phase space region, the soft photon approximation can be used to calculate 
the cross section, provided that 6 S is sufficiently small. The soft singularities are regularized 
by giving the photon a fictitious small mass. In the sum of the virtual and soft photon 



L A parton level Monte Carlo program (in FORTRAN) for pp — ► — ► t^v including 0(a) EW 



corrections is available from [http : //ubhex . physics . buffalo . edu/~baur/wgrad/wgrad . tar . gz 



or by contacting baur@ubhex.physics.buffalo.edu. 
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terms the unphysical photon mass dependence cancels, and the QED-like contributions are 
IR finite. 

The IR finite remainder of the virtual photonic corrections and the pure weak one- 
loop corrections of Fig. [I] can be combined to separately gauge invariant modified weak 
contributions to the W boson production and decay processes. Both the QED-like and the 
modified weak contributions are expressed in terms of form factors, Fq ED and F£ eak , which 
multiply the Born cross section JTH] . The superscript a in the form factors denotes the initial 
state, final state or interference contributions. 

The complete 0(a 3 ) parton level cross section of resonant W production via the Drell- 
Yan mechanism —> £ + u( , ~f) can then be written as follows pi 



+ E [dcr^ F« ED (s,t)+da^ 3 ], (1) 



a— initial , f inal , 
inter f. 



where the Born cross section, da^°\ is of Breit-Wigner form and s and t are the usual 
Mandelstam variables in the parton center of mass frame. The modified weak contributions 
have to be evaluated at s — fTjf] . Explicit expressions for the form factors Fq ED , F£ eak 



are given in Ref . |L8| . The IR finite contribution do r%_^ 3 describes real photon radiation with 

Additional singularities occur when the photon is collinear with one of the charged 
fermions. These collinear singularities are regularized by retaining finite fermion masses. 
Thus, both da an d Fq ED (a = initial, final) contain large mass singular logarithms 
which have to be treated with special care. In the case of final state photon radiation, 
the mass singular logarithms cancel when inclusive observables are considered (KLN theo- 
rem [27|). For exclusive quantities, however, these logarithms can result in large corrections, 
and it may be necessary to perform a resummation of the soft and/or collinear photon emis- 
sion terms. To increase the numerical stability of the inclusive calculation, it is advantageous 
to extract the collinear part from da^z an d perform the cancellation of the mass singular 
logarithms analytically. The reduced 2^3 contribution, i.e. the real photon contribution 
away from the soft and collinear region, can be evaluated numerically using standard Monte 
Carlo techniques. 

For initial state photonic corrections, the mass singular logarithms always survive. These 
logarithmic terms are equivalent to the 1/e singularity encountered in dimensional regular- 
ization (D = 4 — 2e is the number of dimensions) with massless quarks. They are universal to 
all orders in perturbation theory, and can therefore be canceled by universal collinear counter 
terms generated by 'renormalizing' the parton distribution functions (PDF's), in complete 
analogy to gluon emission in QCDQ In addition to the collinear counterterms, finite terms 
can be absorbed into the PDF's, introducing a QED factorization scheme dependence. We 
have carried out our calculation in the QED DIS and QED MS scheme. The extraction of 
the collinear part from d&2^ and the renormalization of the PDF's are described in Sec. IIA 



Alternatively, these logarithmic terms can be retained in the calculation. They would lead to 
large corrections, but then also to large changes in the input PDF's. 
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and IIB, respectively. In Sec. IIC we provide explicit expressions for the 0(a 3 ) cross section 
for W production in hadronic collisions in the QED DIS and QED MS scheme, and study 
the dependence of the 0(a 3 ) cross section on the theoretical cutoff parameters which define 
the soft and collinear regions. 



A. The Extraction of Collinear Singularities from d&2^ 3 

The contribution of real photon emission to the 0(a 3 ) cross section for W production in 
hadronic collisions is given by 

da real = dP 2 ^ 3 ^\M B r\ 2 , (2) 

where <iP2->3 is the product of the three particle phase space element and the flux factor, 

11 d 3 p f d 3 pfd 3 k ,. 
dP ^ =2lj2^ 8 P y fk o + p f - ^ (3) 

and the Bremsstrahlung matrix element M. br is given by 

M BR = i^ r v 7 ^ | 7 _ 1 M ^ u f G^j (1 - 7 5 ) Vf> Vi,Y(l ~ 7sK 

u f 7 „(1 - 7s) v r v v GT(1 - 75) u\ e*{k) . (4) 



s - Myy - 2kq 



In Eq. (f|), = sin 2 8 w , where 9 W is the weak mixing angle, e p denotes the photon polar- 
ization vector, and 



rP = Q (p p f + Y #/2)7 M _ ; r (p p f >+ #772) _ fy+jY-f 1 

kpf kpf kq 



rw _ n r(p P - ^772) n fog-r wvr ^ q p-^Y + g w ft M 

— Hi ; Hv ; ; • V>) 

kpi kpi> kq 

Q a (a = i, i', f, /') denotes the electric charge in units of the proton charge, e. The initial 
and final state currents are separately conserved: k p G^ p = (Qf — Qf — 1) 7 M = and 
k p Gf = {Q l -Q v -l)^ = Q. 

da rea i can be decomposed into soft and hard initial state, final state and interference 
terms: 

da real = J2 + ^2-3)- (6) 

a — initial , f inal , 
inter f. 

Here, da® f t are the soft photon contributions (J5 7 < E cut ) and, as explained before, are 
included in the QED-like form factors Fq ED . 

The initial and final state hard photon contributions, da^^ al and do^Zt ■> contain mass 
singular logarithmic terms, whereas the interference contribution, da^l \ does not. In 
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order to extract the mass singular terms from da 2 _^ 3 
region by requiring that 

cos 6 > 1 — 5a, 



and d<j 2 ^ , we define a collinear 



(7) 



where 9 is the angle between the charged fermion and the emitted photon in the parton 
center of mass frame. da 2 ^ 3 can be decomposed into a finite contribution away from the 
soft and collinear singularity, dd 2 ^ %te , "which will be evaluated numerically, and a collinear 
part do a coll , for which the integration over the singular phase space region can be performed 
analytically, 



i a a i 7 a, finite 

2 ^3 - da colL + da 2 '_i 3 



a = initial, final. 



(8) 



In the following we calculate da^ oll explicitly for both initial and final state photon radiation. 

In the collinear region, dP 2 _^ 3 factorizes into a two particle and a collinear part (see Fig. |2| 
for notation) 



initial state: 

dP 2 ^ 3 (i + i' 
final state: 



f + f+i)->dp 2 ^ 2 (h+t>^f+n 



zd?k 



2(27r) 3 fc° 



-dP 



2^2 



zdzd5i 
16tt 2 



dP, 



2^3 



'l + l' 



f+f+i)->dp 2 _> 2 (i+t'-*h+n 



2 d s k 



2(27T) 3 fc° 



-dp 



2^2 



dzdbj 

16/T 2 



(9) 



(10) 



where we have used 



d A k = 2n(k ) 2 dk°dcos9 



d8ij = 2k u plfd cos 9 ; dk u = -v[dz = ~^p"dz 



(11) 



Using the leading pole approximation, the squared matrix element for initial and final state 
photon emission (see Eq. (f|)) factorizes into the leading-order squared matrix element, 
lA^I 2 , and a collinear factor, c^h-y or c/ 7 ^, provided the parameter 5g is sufficiently 
small: 



initial 1 2 / 
BR I \ 



+ i'^f + f' + i)^Y,\M {0) \ 2 (h + i > 



final 1 2 / 
BR I I 



(0)|2, 



Here 



-i— >h-y 



C f~f^h 



^ 2 1 

7T dj 

,2 « n 2 1 

H 1 Of 



1 + z 2 
z{l-z) 

'1 + z 2 
l-z ' 



i + i' 



2mt 



h + f) c h-*h ■ 



(12) 



2m) 
~J7 



(13) 



and nii (mf) is the mass of the initial (final) state fermion which emits the photon. Com- 
bining Eqs. (|ToD and (|l"2"D , the hard photon contribution in the collinear limit reads (see 
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also Ref. |2§) 
initial state: 



j ~ initial 
d(T coll. 



final state: 



j « final 



1-S S 



l-S s 



dz^da (0) (h + i' ->/ + /' 



dzlda^^Q 2 



(l-z)m^ 
Si'h Sg l\ 



1 111 l 2 o 

1 — z \ mf 2 z 



lQn 2 

2z 



1 - z 



(14) 



l-5s 



dz 



l-Ss 



i + i'^h + f 
1 + z 2 



2z(l-z)(p° h ) 2 S 9 1 



dzlda^—Ql 



1-z 



In 



{X—z)/zm^ 
Sf'h 5q 2 



2^ 
1 - z 



+ (/<-> n 



(15) 



with Si'h-j'h = (p/i + Pi'-j') 2 - I n order to avoid double counting in the soft region, the 
upper limit in the z integration has to be reduced from z = 1 to z = 1 — 5 S in Eqs. ( |T4"D 
and fllBD. 



B. Mass Factorization: QED DIS and MS Scheme 



The mass singular logarithmic terms of Eq. (|14|) can be absorbed by the counter terms 
to the PDF's. In addition to the singular terms, finite 0(a) terms can be absorbed into the 
PDF's. At next-to-leading order (NLO) in QED, the parton distribution functions therefore 
depend on the QED factorization scheme used. In this subsection, we derive the NLO PDF's 
in the QED DIS and the QED WS scheme. 

In order to derive the parton distribution functions at next-to-leading order in QED, 
one must calculate the virtual and real photon contribution to the square of the parton 
electromagnetic current, integrated over the phase space of the final state partons, WL,. 

The contributing Feynman diagrams are shown in Fig. [3]. The tensor WL, is related to the 
structure function F 2 (x, Q 2 ) by (see Fig. ^| for notation) 



F 2 (x,Q' 



x 

47T 



E 



dy 

x y 



I2x 2 

Q 2 ? 



p^ pi' 



w, 



au i 



(16) 



where the sum is taken over all contributing quark flavors, and the qi(y) are the unrenor- 
malized quark distribution functions. Since it involves an additional power of a, we do not 
take into account the photon content of the proton in our calculation. 

In the physical (DIS) scheme the 'renormalized' quark distribution functions are 
defined by requiring that F 2 (x,Q 2 ) is given by the sum of the quark distributions to all 
orders in perturbation theory 



F 2 {x,Q 2 



x 



-.DIS, 



x,Q 2 



+ Q? IS (x,Q' 



(17) 
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where the QED factorization scale has been set equal to Q. The 0(a) structure function 
F 2 (x,Q 2 ) can be obtained from the corresponding 0(a s ) QCD structure function |30| by 
the replacement 



a, 4 



a 



Ql ■ 



71 3 7T 

For massive fermions one finds in the limit Q 2 ^> m 2 (z = x/y): 



(18) 



F 2 (x, Q 2 )=xJ2 Ql { [ f Qi (~ z ) Ul ~z) + l Ql 



Q 2 



21n<5 s + - In ) 5(1 - z 



1 + z 2 , 

+ In 

1 — z 



'q 2 i x 



0(1 -6,-z) 



7 



7T 



\n^5 s + - ln5 s + - + — 



8(1 -z) + 



1 1 



2 1-z 



3z 



0(1 -5,-z) 



• (19) 



Using Eq. ( |l9|) it is then straightforward to calculate qf IS (x,Q 2 ) in terms of the unrenor- 
malized quark distribution functions q%(x) (see below). The relation between F 2 (x, Q 2 ) and 



the quark distribution functions in the MS scheme |3T[ is given by 
F 2 (x, Q 2 )= X J2 Ql [<lF& Q 2 ) + Q 2 

i 

'i dy 



x y 



1 a I 




-cA 


i) 


J 7T \ 





(20) 



with 32 



Ci(z) = -Q\ 



In o s mo s 1 



+ 



1 + z 2 1-z 3 1 

In h 2^ + 3 



8(1 -z) 

1-8 S - z) 



(21) 



l-^ z 21-z 

Ci(z) represents the finite part of the QED 0(a) corrections to deep inelastic scattering after 
removing the singularities according to the MS prescription. The MS scheme is defined in 
the framework of dimensional regularization but Eq. ( p0[ ) can also be used for its definition. 
To obtain the 'renormalized' quark distribution functions, q^ (x,Q 2 ), when finite quark 
masses are used as regulators, we make use of the relation 



MSi 



dy 



x, Q 2 ) + I ^ qf S (V, Q' A ) -4-)= q? 1 *^ Q 



a x 



DIS I 



x y 



(22) 



which follows from Eqs. fll7|) and (|20|) . 

The final expression for the scheme dependent 'renormalized' quark distribution function 
in NLO QED is 
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qi(x,Q 2 ) = qi(x) 



a 



1 + - QUI - \nS s - \n z 5 S + [\n8 s + - ) In 



71 



Q 2 



+ 



i- s <> dz fx\ a 2 ( 1 + z 2 (Q 



x Z 



zj 2tt 



111 



1 — z \m 2 (1 — z) 2 



1 + z 2 
1 - z 



+ \ FV f c \ , (23) 



with 



and 



27T 2 

f v+s = 9 + — + 3\n5 s -2\n 2 5 s , 



, 1 + z 2 ^ fl-z\ 3 1 
/c = In 



l-z 



2 1 - z 



2z + 3 . 



(24) 



(25) 



The QED D/S 1 (MS") scheme corresponds to = 1 (X^c = 0). 



C. The Cross Section for pp — » W(7) — ► ^(7) 



The differential cross section for pp — > ^(7) — » ^(7) is obtained by convoluting the 
parton cross section of Eq. (|T]) with qi(x), and subsequently replacing the unrenormalized 
quark distribution functions by qi(x,Q 2 ), using Eq. (|23|). The initial state QED-like contri- 



bution da^ Fqed 0,1 and the collinear part da l ™[l ta \ including the effect of mass factorization, 
can be grouped into a single 2^2 contribution: 



d(J initial = Y^J dxidx2 [ qi ( Xl ,Q 2 ) q it {x 2 ,Q 2 ) dff® + (1 



2)] 



l,V 



a 



71 



ln5 s + 



In 



, x 3 «* 
\n5 s H 1 

2 24 



Tl 2 I 

+ — - 2 + In 2 5, + -Xpc f vA 
4 



+ ^2 J dx x dx 2 I J 

i.i' (. x 



da® 



a 
27 



l ~ s * dz 

1 + Z 2 ' 



Q 2 ft 



1 



In 



z 



Ail 

Q 2 



X-2 



Q 2 ) frOci, Q 2 ) + Q 2 gi(ar ls Q 2 ) fc, Q 2 



) 2 6* 



+ l-z- Xpcfc + (1 <-> 2) 



(26) 



As expected, the mass singular logarithms cancel completely. x\ and 22 in Eq. (26) are the 
momentum fractions of the parent hadrons carried by the partons. 

In order to treat the 0(a) initial state QED-like corrections to W production in hadronic 
collisions in a consistent way, QED corrections should be incorporated in the global fitting 
of the PDF's using the same factorization scheme which has been employed to calculate the 
cross section. Current fits [33] to the PDF's do not include QED corrections. A study of the 



effect of QED corrections on the evolution of the parton distribution functions indicates [19 
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that the modification of the PDF's is small. We have not attempted to include QED correc- 
tions to the PDF evolution in the calculation presented here. The missing QED corrections 
to the PDF introduce an uncertainty which, however, is likely to be smaller than the present 
uncertainties on the parton distribution functions. 

The squared matrix elements for different QED factorization schemes differ by the finite 
0(a) terms which are absorbed into the PDF's in addition to the singular terms. In the QED 
DIS scheme, the contribution of the QED-like initial state corrections to the cross section 
is about 8% smaller than in the QED MS scheme. The factorization scheme dependence is 
expected to be reduced when the O(o) QED corrections to the PDF are included. In the 
following, for the numerical evaluation of da™^ 0,1 we use the QED MS scheme. 

The final state 2-^2 contribution can be obtained directly from the form factor Fq^ 1 
of Ref. |H| with Qf =u = and Qp=e = —1: 

d4Z2 = E/ dx ^ M^Q 2 ) Mx2,Q 2 ) da^ + (1 <-> 2)] 



In sufficiently inclusive observables the mass singular logarithmic terms cancel in the sum 

r j final j j final 
Of d<7 2 ^2 an d «02-*3 • 

The approximation used so far in modeling the EW radiative corrections to W boson 



production at the Tevatron |13| , |14|| ignores all weak, interference and initial state photonic 
corrections, and differs from our calculation in the treatment of the final state virtual and 
soft photon contribution. At the parton level, the difference between Eq. (p7j) and the 2-^2 



contribution to the differential cross section in the approximate calculation, is given by 

In Sec. IIIA we shall demonstrate that the difference has a non-negligible effect on the shape 
of the transverse mass distribution. 

Experimentally, photons which are collinear with muons can be identified without prob- 
lems: photons deposit energy in the electromagnetic calorimeter, whereas muons are iden- 
tified by hits in the muon chambers. For muons in the final state, therefore, one has to 
retain full information on the particle momentum four vectors. In the electron case, on the 
other hand, the finite resolution of detectors makes it difficult to discriminate between elec- 
trons and photons with a small opening angle, and the electron and photon four-momentum 
vectors are recombined to an effective electron four-momentum vector if their separation 
Ai? e7 in the azimuthal angle-pseudo-rapidity plane is smaller than a critical value R c . If 
the lepton and photon four-momentum vectors are not resolved in the collinear region, the 
collinear singularities from the hard photon contribution can be extracted as described in 
Sec. IIA, and the integration over the momentum fraction z in Eq. ( |1~5| ) can be performed 
analytically. The parameter Sg has to be chosen sufficiently small to ensure that AR ei < R c 
over the entire region where the analytic integration is carried out. For small R c one finds 
that 5g has to be less than 

d2 

Aiiia.N ^ (29) 



2coslr(77 max (e)) 
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where r] max (e) is the maximum allowed pseudo-rapidity of the electron. 

The procedure described above is part of the electron identification process used by the 
D0 collaboration pi. The CDF collaboration uses a slightly different method where the 
electron and photon four-momentum vectors are combined if both particles traverse the 
same calorimeter cell. This modifies the expression for 5™ ax . 

Once the integration over z has been performed analytically, the collinear contribution 
dc r ioU, al an d the QED-like contribution da^Fg^ can be combined to cancel the mass sin- 
gular logarithms explicitly, and one obtains for the final state 2 —>■ 2 contribution 



^^final 



'2^2 



J dx\dx 2 \qi(xi,Q 2 ) 



2)] 



a. 
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7T 



x - { -In (J. + - - ln<L + - In -f 



7T 



4 



4 



'5, 



(30) 



As described in Ref. ||18|1 , a part of the photonic interference terms together with the IR 
finite parts of the box diagrams in Fig. |1] can be absorbed into the modified weak contribu- 
tions p™^ l 'f mal introduced earlier in this Section. The 2^2 interference contribution is 
then given by 



d^nterf. _ J dx x dx 2 qi(xi,Q 2 ) q i ,(x 2 ,Q 2 ) <ia (0) 



x -f3 int .{s,t,u) In 



M^-s5 s y + M^r 2 Wi 



2) 



(31) 



where Y w is the W width, and 



a. 



7T 



Q v \n{ +2 



(32) 



For s = Myy, the 2^2 interference contribution is completely cancelled by the hard photon 
contribution da^l when the total inclusive cross section is calculated. Note that the 2-^2 
interference contribution exhibits only soft singularities. 

The complete 0{p?) cross section for pp — > W{^) — > tv^) can now be expressed as 

^(o+D = Y,J dx,dx 2 MxuQ 2 ) qA*2,Q 2 ) da® + (1 <-> 2)] 



initial 



x [1 + 2 Tie (F weak 

initi 
r 2^3 



final 



weak 



){M 2 W )\ 



+ daiZf 1 + daf*i alJinite + da^ e 2 rf + daf^ + da fmal 



(33) 



with 



da final = d4Z 2 + dalZf 



(34) 



if the integration over z is performed numerically, and 

da fmal = da f 2 ™ 2 + dag!g Jinite 



(35) 
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if the z-integration is done analytically. Here, da^l?™^ (a = initial, final) are the reduced 
2^3 contributions away from the soft and collinear region. The hard bremsstrahlung 
contribution has been compared numerically with the pp — > £wy cross section of Ref. [34 |. 
The two calculations agree to better than 1%. 

The end result of the calculation consists of two sets of weighted events corresponding 
to the 2 — > 2 and 2-^3 contributions. Each set depends on the parameters 5 S and Sq. 
The sum of the two contributions, however, must be independent of 5 S and Sg, as long 
as the two parameters are taken small enough so that the soft photon and the leading 
pole approximation are valid. In Figs. [| - |6] we show the different contributions to the 
pp —>■ £ + i>('-f) cross section at -v/s = 1.8 TeV as a function of the two parameters. To 
compute the cross section, we use here and in all subsequent figures the MRSA set of parton 



distribution functions [ 35f| , and take the QCD renormalization scale uqcd and the QED and 



QCD factorization scales, Mqed and Mqcd, to be Uq CD = Q 2 = Mq ED = Mq CD = M^- 
The detector acceptance is simulated by imposing the following transverse momentum (pt) 
and pseudo-rapidity (n) cuts: 

p T (£) > 25 GeV, \n(£)\ < 1.2, £ = e, u, (36) 



i> T > 25 GeV. (37) 

These cuts approximately model the acceptance cuts used by the CDF and D0 collaborations 
in their W mass analyses []T5| , p!6| . Uncertainties in the energy and momentum measurements 



of the charged leptons in the detector are simulated in the calculation by Gaussian smearing 
of the particle four-momentum vector with standard deviation a which depends on the 
particle type and the detector. The numerical results presented here were calculated using 



a values based on the specifications for the upgraded Run II D0 detector p6| . The results 
obtained using the target specifications for the CDF II detector [^] are similar. The SM 
parameters used in our numerical simulations are M w = 80.3 GeV, M z = 91.187 GeV, 
a = a (o) = 1/137.036, G M = 1.166 x 10~ 5 GeV~ 2 , T w = 2.1 GeV, and m top = 175 GeV. 
These values are consistent with recent measurements at LEP, SLC and the Tevatron M. 

Figure |] displays the QED-like initial state (ISR) corrections to the cross section as a 
function of 5 S (Fig. f|a) and 6$ (Fig. [|b). In order to exhibit the independence of the cross 
section from the parameters S s and Sg more clearly, we have not included the Born cross 
section in the 2-^2 contribution here as well as in Figs. § and ||]. The QED-like ISR 
corrections to the cross section for electron and muon final states are virtually identical. 
While the separate 2^2 and 2^3 0(a) contributions vary strongly with 6 S and 5q, 
the sum is independent of the two parameters within the accuracy of the Monte Carlo 
integration. 

In Fig. we show the QED-like final state (FSR) corrections to the pp — > £ + i , ('-f) cross 
section as a function of 5 S for electron and muon final states. Radiation of photons collinear 
with one of the leptons gives rise to terms proportional to [ln(s/m 2 ) — 2] \n(S s ) (see Eq. (|2"7D) 
in both the 2 — > 2 and 2 — > 3 contributions. As demonstrated in Fig. |5], these terms cancel 
and the total cross section is independent of 5 S . Due to the smaller mass of the electron, the 
variation of the 2 — > 2 and 2 — > 3 contributions with 5 S is more pronounced in the electron 
case. 
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In Figs. f| and [5], we have not taken into account realistic lepton identification require- 
ments, i.e. we have assumed that photons and leptons with arbitrary small opening angles 
can be discriminated. In a more realistic simulation, in addition to the lepton p?, f T and 
pseudo-rapidity cuts, one imposes requirements on the separation of the charged lepton and 
the photon. These requirements differ slightly for the CDF and D0 detectors. In the fol- 
lowing we adopt lepton identification criteria which are motivated by the D0 W mass [16| 



and W cross section PS| analyses; the numerical results obtained using the requirements 
imposed in the CDF W mass analysis |15j are similar. In order to study their impact on 



the size of the EW radiative corrections, we will perform simulations both with and without 
the lepton identification requirements taken into account. 

We shall use the following lepton identification requirements. For electrons, we require 
that the electron and photon momentum four-vectors are combined into an effective electron 
momentum four-vector if AR n < 0.2. For 0.2 < AR n < 0.4 events are rejected if E n > 
0.15 E e . Here E 1 (E e ) is the energy of the photon (electron) in the laboratory frame. For 
events with 0.2 < Ai? e7 < 0.3 and i? 7 < 0.15 E e , the electron and photon momentum four- 
vectors are again combined. Muons are identified by hits in the muon chambers and the 
requirement that the associated track is consistent with a minimum ionizing particle. This 
limits the photon energy for small muon - photon opening angles. For muons, we therefore 
require the energy of the photon to be E 1 < 2 GeV for A_R M7 < 0.2, and E 1 < 6 GeV for 
0.2 < Ai2 M7 < 0.6. For future reference, we summarize the lepton identification requirements 
in Table |. 

As noted before, when the electron and photon momentum four-vectors are combined, it 
is possible to analytically cancel the mass singular terms in the QED-like 2 — > 2 final state 
corrections. In this case, the QED-like final state 2 — > 2 and 2 — > 3 contributions depend on 
the collinear cutoff parameter 8g. Figure |5] displays the QED-like final state contribution to 
the pp — » e + z/(7) cross section as a function of 5$ when the electron identification require- 
ments described above are taken into account. While the 2 — > 2 and 2^3 contributions 
both exhibit a considerable dependence on Sg, their sum is independent of the parameter. 

Similar to the QED-like initial and final state corrections, one can show that the sum of 
the 2^2 and 2^3 contributions of the QED-like initial - final state interference terms 
is independent of 5 S . The interference terms are typically of the same size as the initial 
state corrections. The modified weak contributions to the 0(a 3 ) cross section are trivially 
independent of 5 S and 5g. In the following, these parameters will be fixed to 5 S = 10 -2 and 
5 e = 10- 3 . 

As stated before, we take the QCD renormalization scale fiQCD and the QED and QCD 
factorization scales, Mq ED and Mq C d, to be equal, Hqcd = Mq ED = Mq C d = Q- The 
missing QED corrections to the PDF's create a dependence of the O(o) initial state cor- 
rections on the scale Q which is stronger than that of the lowest order calculation. On the 
other hand, final state and initial - final state interference terms depend on Q only through 
the PDF's. 

III. PHENOMENOLOGICAL RESULTS 

We shall now discuss the phenomenological implications of the 0(a) electroweak correc- 
tions to W production at the Tevatron (pp collisions at y/s = 1.8 TeV). We first discuss the 
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impact of electroweak corrections on observables used to measure Mw' the transverse mass 
distribution, the Pt{£) distribution, and the W to Z transverse mass ratio. We then consider 
the W production cross section, the W to Z cross section ratio and the charge asymmetry of 
leptons in the W decay. Unless stated otherwise, we take into account the cuts of Eqs. ( p6|) 
and (p7| ) and effects from energy and momentum measurement uncertainties in the detector. 
We state explicitly when the lepton identification requirements listed in Table | are included. 



A. Electroweak Corrections to the Mr and pr(() Distributions, and the W to Z 

Transverse Mass Ratio 



Since the detectors at the Fermilab Tevatron collider cannot directly detect the neutrinos 
produced in the leptonic W boson decays, W — > iv, and cannot measure the longitudinal 
component of the recoil momentum, there is insufficient information to reconstruct the 
invariant mass of the W boson. Instead, the transverse mass distribution of the final state 



lepton pair, or the transverse momentum distribution of the charged lepton are used [16| to 
extract Mw- The transverse mass is defined by 



M T = pp T (£)p T (v)(l - cos0^) , (38) 

where Pt{P) and priy) are the transverse momentum of the lepton and the neutrino, and 
(f) is the angle between the charged lepton and the neutrino in the transverse plane. The 
neutrino transverse momentum is identified with the missing transverse momentum, p" T , in 
the event. Recently, it has been pointed out that the ratio of W to Z observables can 
also be used to measure the W mass [[23 . This method has been applied to the W to Z 



transverse mass ratio by the D0 collaboration ||24j| . The advantages and disadvantages of 
the observables used to extract M w are discussed in Ref. |J. 

The 0(a 3 ) Mt distribution for e + z/(7) (solid) and fi + u('~f) (dots) production is shown 
in Fig. 0a together with the lowest order predictions (dashed and dot-dashed curves). In 
Fig. ^|b, we show the O (a 3 ) and Born Pt{£) spectrum. The flavor specific lepton identification 
requirements of Table | are not taken into account here. Electroweak corrections decrease 
the cross section at the peak of the M T (p T (£)) distribution by about 12% (17%) in the 
electron, and by about 6% (7%) in the muon case. Photon radiation from the charged 
lepton lowers the Iv invariant mass. Events from the Jacobian peak regions in the Mr and 
Pt (£) distributions therefore are shifted on average to lower values of the transverse mass and 
transverse momentum. Due to the ln(s/mf ) term, the effect of the corrections is larger in the 
electron case. The Jacobian peak of the pt distribution is broader and less pronounced in 
the muon case, due to the energy and momentum resolution which is significantly worse for 
muons than for electrons. In the Mr distribution, the effect of the j) T resolution dominates, 
and the difference in the distribution between electrons and muons is small. The Mr and 
Pt{£) distributions for £~u('~f) production are identical to those for the £ + z/(7) channel in pp 
collisions. In the remainder of this subsection we therefore only consider the i^vi^j) final 
state. 

The various individual contributions to the EW O{o) corrections on the Mr distribution 
are shown in Fig. The initial state QED-like contribution uniformly increases the cross 
section by about 1% for electron (Fig. ||a) and muon (Fig. ||b) final states. It is largely 
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canceled by the modified weak initial state contribution. The interference contribution is 
very small. It decreases the cross section by about 0.01% for transverse masses below My/, 
and by up to 0.5% for M? > My/. The final state QED-like contribution significantly 
changes the shape of the transverse mass distribution and reaches its maximum effect in the 
region of the Jacobian peak, M? ~ My/. As for the initial state, the modified weak final 
state contribution reduces the cross section by about 1%, and has no effect on the shape of 
the transverse mass distribution. For Mt > 125 GeV, the QED-like final state corrections 
uniformly reduce the differential cross section by about 5% in the electron case, and by about 
2% in the muon case. Without taking the lepton identification requirements of Table [I] into 
account, the full O{o) electroweak radiative corrections to the Mt distribution are very well 
approximated by the sum of the QED-like and modified weak final state corrections. 

It should be noted that the differential cross section ratio shown in Fig. [| becomes ill 
defined in the threshold region Mt ~ p'r' t {€) + ftp , where p^*(£) and are the charged 
lepton pt and the missing transverse momentum threshold. For Mt < p T ut {€) + , the 
Born cross section vanishes, and the cross section ratio is infinite. The 0(a 3 ) cross section 
is small, but non-zero, in this region. The largest contribution to the cross section for 
Mt < Pt" 1 ^) +^T* t originates from initial state radiation configurations, where the lepton 
and the neutrino have a small relative opening angle and are balanced by a high pt photon 
in the opposite hemisphere. Close to the threshold, Mt ~ p^if) +p^ t , large logarithmic 
corrections are present, and for an accurate prediction in this region those corrections need 
to be resummed. The results of Fig. [8] in this region should therefore be interpreted with 
caution. 

The ratio of the full 0(a 3 ) and the Born cross section as a function of the transverse 
mass is shown in Fig. || The solid (dashed) lines show the cross section ratio without 
(with) the lepton identification requirements included. Recombining the electron and pho- 
ton four-momentum vectors for Ai? e7 < 0.2 eliminates the mass singular logarithmic terms 
and strongly reduces the size of the QED-like final state corrections (see Fig. |S]a). These 
corrections are now of the same size as the initial state QED-like and the modified weak 
corrections. However, with the total 0(a) EW corrections varying between 1% and 2%, 
the shape change of the Mt distribution caused by the final state corrections is still signif- 
icant. For muon final states (see Fig. ||b), the cut on the energy of the photon reduces the 
hard photon part of the 0(a 3 ) ^(7) cross section. In this case, the mass singular terms 
survive and the corrections become larger over the entire range of Mt considered. Before 
lepton identification requirements are taken into account, the change in the shape of the Mt 
distribution due to the QED-like final state corrections is more pronounced in the electron 
channel. Once these requirements are included, the shape change is stronger in the muon 
case. 

Results qualitatively similar to those shown in Figs. |8] and ^ are obtained for the trans- 
verse momentum distribution of the charged lepton. 

As we have seen, final state bremsstrahlung has a non- negligible effect on the Mt and 
Pt{£) distribution in the Jacobian peak region. As is well known, electroweak corrections 
must be included when the W boson mass is extracted from data, otherwise the mass is 
shifted to a lower value. In the approximate treatment of the electroweak corrections used 
so far by the Tevatron experiments, only final state QED corrections are taken into account; 
initial state, interference, and weak correction terms are ignored. Furthermore, the effect of 
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the final state soft and virtual photonic corrections is estimated from the inclusive 0(a 2 ) 
W — > tvid) width |J8| and the hard photon bremsstrahlung contribution [13,14]. When 
detector effects are included, the approximate calculation leads to a shift of about —50 MeV 
in the electron case, and approximately —160 MeV in the muon case |15|J16[ . 



Initial state and interference contributions do not change the shape of the Mt distribution 
significantly (see Fig. |8|) and therefore have little effect on the extracted mass. However, 
correctly incorporating the final state virtual and soft photonic corrections results in a non- 
negligible modification of the shape of the transverse mass distribution. This is demonstrated 
in Fig. [IT], which shows the ratio of the Mt distribution obtained with the QED-like final 
state correction part of our calculation to the one obtained using the approximation of 
Refs. 11 and @. The depend ence of the ratio on Mt is described by Eq. (|28|). For Mt < 
Mw, most events originate from the region s « M w , due to the Breit-Wi gner resonance. 
Consequently, there is little dependence of the cross section ratio on Mt in this region. For 
Mt > M\y, the steeply falling cross section in the tail of the Breit-Wigner resonance favors 
events with s ~ M T . In this region the term proportional to \n(s/M w ) in Eq. (^) causes a 
change in the shape of the transverse mass distribution. Aa^ mal also contains a term which 
is proportional to \ia(mj/M w ) (see Eq. (|28])). The shape change in the Mt distribution 
thus is more pronounced in the electron case. Lepton identification requirements have a 



small effect on the cross section ratio (see Fig. 10). Note that the approximate NLO cross 
section, and thus the cross section ratio shown in Fig. [l(| does depend explicitly on the 
cutoff S s whereas the 0(a 3 ) cross section resulting from our calculation does not. While the 
dependence on 5 S is very small for Mt < Mw, it is quite pronounced for transverse masses 
above Mw- 

The difference in the line shape of the Mt distribution between the complete 0(a 3 ) 
calculation and the approximation used so far occurs in a region which is important for both 
the determination of the W mass, and the direct measurement of the W width. The precision 
which can be achieved in a measurement of Mw using the transverse mass distribution 
strongly depends on how steeply the Mt distribution falls in the region Mt ~ Mw (see 
Fig. [?]). In the region of large transverse masses, My > 100 ... 110 GeV, the shape of the 
Mt distribution is sensitive to the W width |[39| . Any change in the theoretical prediction 
of the line shape thus directly influences the W mass and width measurements. From a 
maximum likelihood analysis similar to that carried out in Ref. [^TJ for Z production, the 
shift in the measured W mass due to the correct treatment of the final state virtual and 
soft photonic corrections is found to be AMw ~ O(10 MeV). This shift is much smaller 
than the present uncertainty for Mw from hadron collider experiments |15],|T6|. However, for 



future precision experiments, a difference of 0(10 MeV) in the extracted value of Mw can 
no longer be ignored, and the complete O (a 3 ) calculation should be used. 

At high luminosities, the transverse mass ratio of W to Z bosons offers advantages in 



determining the W mass p3| , |24| over the Mt and pr(^) distributions. The transverse mass 
ratio of W and Z bosons is defined as 



Rm t {X Mt ) 



Xm 7 



A Z (X Z 



M T 



Xm 7 



(39) 



where Ay (V = W, Z) is the differential cross section 
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Mx v Mt ) = ^- do) 



with respect to the scaled transverse mass, 



X v - M¥ (41) 



The transverse mass of the lepton pair in Z boson events is defined in complete analogy to 
Eq. ©: 

M| = V / 2j9 T (£ + )pt(^-)(1-cos0) , (42) 

where is the angle between the two charged leptons in the transverse plane. 

The ratio of the 0(a 3 ) and the Born W to Z transverse mass ratio is shown in Fig. 11. 
To calculate the 0(a) electroweak corrections to Z boson production, we use the results of 



Ref. [21]. Note that purely weak corrections are not included in this calculation. Identical 
charged lepton and rapidity cuts are used for W and Z production. In the Z boson 
case, photon exchange and jZ interference effects are included and an additional cut on the 
di-lepton invariant mass of 75 GeV < m(l + £r) < 105 GeV has been imposed. In pp — > ivi^f) 
only one of the two leptons can emit a photon, whereas both leptons in pp —>■ £ + £~ (7) can 
radiate. The 0(a) corrections are thus significantly larger in the Z case. As a result, the W 
to Z transverse mass ratio is more strongly affected by electroweak radiative corrections than 
the Mjf distribution. Without the lepton identification requirements, the O(o) corrections 
increase Rm t by about 30% (10%) at the location of the Jacobian peak (X Mt = 1) for 
electrons (muons). For X Mt < 0.9, the electroweak corrections reduce the transverse mass 
ratio by 6 — 10% in the electron case, and by 4 — 6% in the muon case. 

When lepton identification criteria are taken into account, the merging of the electron 
and photon momentum four- vectors for small e — 7 opening angles again strongly reduces 



the size of the O(o) corrections (see Fig. Ob). In the region of the Jacobian peak, the 



corrections are reduced to ~ 4%, and for X Mt < 0.95 to about 2% in magnitude. For muon 
final states, the lepton identification requirements reduce the hard photon part of the 0(a 3 ) 
cross sections below the Jacobian peak, but have little effect in the peak region where only 
few events contain hard photons (see Fig. |S|b). The reduction is more pronounced for the 
pp — > n + n~(~f) than for the pp — > ^(7) cross section. Consequently, the 0(a) corrections 
increase Rm t below the Jacobian peak and leave it almost unchanged in the peak region. 



B. Electroweak Corrections to the W Boson Cross Section and the W to Z Cross 

Section Ratio 



In the past, the measurement of the W and Z boson cross sections has provided a test of 
perturbative QCD [[25|,[|(],[|l|] . With the large data set accumulated in the 1994-95 Tevatron 
collider run, the uncertainty associated with the integrated luminosity (~ 3.6% ETJ) became 
a limiting factor in this measurement. This suggests to use the measured W and Z boson 
cross sections to determine the integrated luminosity in Run II f41.42]|. The cross section 
ratio 
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_ a(pp ^W^ luX) 

Rw/Z ~ a( P p ^z^ m-x) ' (b) 



together with the theoretical prediction for the ratio of the total W and Z production 
cross sections, (Jwl&z — 3.36 ± 0.02 [fl3|| , the LEP measurement of the branching ratio 
B(Z — > £ + £~) and the SM prediction for the W —>■ iv decay width, can be used for an 
indirect determination of IV p5lj44|| . F° r integrated luminosities smaller than about 20 fb~ , 
the W width measurement from Rw/z is expected to yield better results than the direct 
determination from the Mt distribution |5j . 

The size of the 0(a) electroweak corrections to the total pp — > IvX cross section and to 
Rw/z is sensitive to the acceptance cuts and whether lepton identification requirements are 
taken into account or not. In Table ||, we list the electroweak .fT-factor, 

REW = a°W(pp^W^evX) 

0-Born(pp _, W _^ ^ ' ^ > 

and the correction factor for Rw/z, 

AT" = ^ ■ («) 

^W/Z 

for the acceptance cuts listed in Eqs. ( |36"D and ( |3"TD with and without taking the lepton 
identification requirements of Table | into account. As before, we include photon exchange 
and 7Z interference effects, and impose a cut on the di-lepton invariant mass of 75 GeV < 
m(l + r) < 105 GeV in the calculation of the 0(a 3 ) Z boson cross section entering 



Rw/z- It should be noted that the missing purely weak corrections in the calculation of 
Ref. [^TJ introduce an uncertainty of 0(a/ir) in K^ w which could be significant. 



From the results listed in Table |T| we see that the 0(a) electroweak corrections decrease 
the W cross section and increase Rw/z by several per cent for the cuts imposed. As for the 
differential cross section, the 0(a) corrections are larger in the electron case when lepton 
identification requirements are not included. When lepton identification requirements are 
included, the corrections are reduced in the electron case and enhanced in the muon case. If 
no acceptance cuts and no lepton identification requirements are taken into account, all mass 
singular terms cancel in the total cross section, and the size of the electroweak corrections 
is reduced to about —0.2%. 

The size of the 0(a) electroweak corrections should be compared with that of the 0(a s ) 
and 0(a 2 s ) QCD corrections. NLO QCD corrections are known (see e.g. Ref. j|5| and 
references therein) to enhance the W production cross section by about 15 — 20% and 
are thus significantly larger than the 0(a) EW corrections. In fact, the size of the 0(a) 
electroweak corrections to the W cross section when cuts are imposed is about equal to that 
of the next-to- next-to-leading order (NNLO) QCD corrections |[46|| . On the other hand, the 



O(a) electroweak corrections to Rw/z are in some cases considerably larger than the NLO 
QCD corrections. Since the QCD corrections to W and Z production are very similar, they 
cancel almost perfectly in the W to Z cross section ratio; the 0(a B ) corrections to Rw/z are 



of 0(1%) or less, depending on the set of parton distribution functions used |jl6"|| . In contrast, 



the electroweak corrections do in general not cancel in Rw/z- As noted before, in Z — > i + i 
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both leptons can emit photons, whereas only the charged lepton radiates in W — > £v decays. 
Since final state photonic corrections are the dominating contribution to the 0{a) EW 
corrections, the 0(a) corrections to the W and Z cross sections are quite different, and 
thus do not cancel in Rw/z- F° r example, when lepton identification requirements are taken 
into account, the 0(a) EW corrections in the electron (muon) case increase Rw/z by 0.2% 
(6.5%). Note that, unlike the electroweak corrections, QCD corrections are only slightly 
modified by cuts and lepton identification requirements. 

From the results shown in Table || we conclude that it will be necessary to correct for 
higher order electroweak effects if one wishes to measure the W cross section and Rw/z with 
an accuracy of (9(1%) or better. 



C. Electroweak Corrections to the Charge Asymmetry of Leptons in W decays 

Uncertainties in the parton distribution functions are a major contribution to the sys- 



tematic error of the W mass extracted in hadron collider experiments |[3,|I|]. Measurement 
of the charge asymmetry of leptons in W decays ^\ , 



da + /dy(£)-da-/dy(£) 
A{M) = da+/dy(l) + da-/dy(l) ' (46) 

where y(£) is the lepton rapidity and 

a ± = a{jrp^l ± vX), (47) 



provides strong constraints on the ratio of d and u quark distributions . These constraints 
considerably reduce the uncertainty originating from the parton distribution functions in 
the W mass measurement JT5|JT6[] . It is thus important to know how electroweak radiative 



corrections affect A(y(£)). The 0(a 3 ) asymmetry as a function of the lepton rapidity for 
eiy(j) (dashed line) and ^(j) (dotted line) production is shown in Fig. [12| together with the 
lowest order prediction (solid line). Except for the pseudo-rapidity cut on the charged lepton, 
we impose the cuts listed in Eqs. (^6]) and fl3T|) in this subsection. Since A(—y) = —A(y), 
the asymmetry is only displayed for y(£) > 0. The flavor specific lepton identification 
requirements of Table | are not taken into account in Fig. The asymmetry in the Born 
approximation for electron and muon final states is virtually indistinguishable for the cuts 
and the energy and momentum resolutions we use. Electroweak corrections are seen to only 
slightly affect the charge asymmetry. 

In order to display the effect of EW radiative corrections on A(y(£)) more clearly, we 
show the difference between the 0(a 3 ) and the Born asymmetry in Fig. [13|. Without taking 
the lepton identification requirements of Table | into account, the difference of the 0(a 3 ) 
and the Born asymmetry is positive and gradually increases with y(£) from zero at y{£) = 
to about 0.01 for electrons, and approximately 0.005 for muons, at y{£) = 2.5 (see Fig. |13|a). 
Figure [T3|b displays the difference between the 0(a 3 ) and the Born charge asymmetry when 
the lepton identification requirements are included in the simulation. Due to the recombina- 
tion of the electron and photon momentum four vectors for small lepton - photon opening 
angles, the size of the electroweak corrections to the charge asymmetry in the electron case 
is drastically reduced when these requirements are taken into account. For y(e) < 1.5, the 
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O(a) corrections reduce A(y(e)) by a very small amount. In the forward rapidity region, 
1.7 < y(e) < 2.5, the difference of the 0(a 3 ) and the Born charge asymmetry is positive 
and slowly increases, reaching approximately 0.0025 at y(e) = 2.5. For comparison, the 
statistical error in A(y(e)) in this region expected for Run II (assuming JCdt = 2 fb _1 ) 
is SA(y = 2.5) ~ 0.005 [^,0. The variation of the charge asymmetry due to the uncer- 
tainties of the present parton distribution functions is about 0.015 {49] in the same region. 
The magnitude of the 0(a) electroweak corrections for muons and electrons when lepton 
identification requirements are included is similar. However, in the muon case, the EW 
corrections enhance (reduce) the 0(a 3 ) asymmetry at small (large) rapidities. 



IV. CONCLUSIONS 

The mass of the W boson is one of the fundamental parameters of the SM and a precise 
measurement of M\y is an important objective for current experiments at LEP2 and future 
experiments at the Tevatron. A precise measurement of the W mass helps to constrain 
the Higgs boson mass from radiative corrections. It will also provide restrictions on the 
parameters of the MSSM. In order to perform such a measurement at a hadron collider, 
it is crucial to fully control higher order QCD and electroweak corrections. In this paper 
we have presented a calculation of the electroweak corrections to W production in hadronic 
collisions which is based on the full set of contributing 0(a 3 ) Feynman diagrams. 

The 0(a) electroweak corrections can be arranged into separately gauge invariant QED- 
like contributions corresponding to initial state, final state and interference corrections, and 
gauge invariant modified weak contributions to the W production and decay processes. Due 
to mass singular logarithmic terms associated with final state photon radiation in the limit 
where the photon is collinear with one of the leptons, final state radiation effects domi- 
nate. Initial state corrections were found to be small after appropriately factorizing the 
corresponding collinear singularities into the parton distribution functions. However, cur- 
rently no parton distribution functions which include QED corrections are available. With 
the factorization scheme used in this paper, the effect of the QED corrections on the PDF 
is expected to be small. We find that the part of the initial state corrections included in 
our calculation is uniform over the entire range of the iv transverse mass and the lepton 
Pt range, and increases the differential cross section by about 1%. Likewise, the modified 
weak corrections are uniform, but decrease the cross section by a similar amount. In con- 
trast, the final state QED-like corrections modify the shape of the and p? distributions 
substantially. 

Without including the lepton identification requirements imposed by experiments, the 
effect of the electroweak corrections is larger in the electron channel than in the muon 
channel. When these requirements are taken into account, the mass singular logarithmic 
terms are eliminated in the electron case because the electron and photon momentum four 
vectors are combined for small opening angles where it is difficult to resolve the two particles. 
Initial state QED-like, final state QED-like, modified weak and interference contributions 
are then all of similar size. On the other hand, in order to experimentally identify muons, 
the energy of the photon is required to be smaller than a critical value if the \i — 7 separation 
is small, and mass singular terms survive. Removing energetic photons thus enhances the 
effect of the O(a) corrections, and the effect of the electroweak corrections in the muon case 
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is larger than in the electron case once lepton identification requirements are included. 

Electroweak radiative corrections have a significant impact on the W mass extracted from 
experiment. The main effect is caused by final state photon radiation and is corrected for in 
the W mass analyses of the Tevatron experiments ||T5| , p!6| . However, in the calculation used 
by CDF and D0 jTSyTJ], the effect of the final state soft and virtual photonic corrections 
is estimated indirectly from the inclusive 0(a 2 ) W — > ivici) width and the hard photon 
bremsstrahlung contribution. Initial state, interference, and weak contributions to the 0(a) 
corrections are ignored altogether. The correct treatment of the final state soft and virtual 
photonic corrections significantly changes the slope of the transverse mass distribution in 
the region Mt > Mw This changes the W mass extracted from the transverse mass 
distribution by (9(10 MeV), and might also have a non-negligible effect on the W width 
measured from the tail of the Mt distribution. More detailed numerical simulations are 
needed to quantitatively assess this effect. Initial state, and initial - final state interference 
corrections, have only a small effect on the Mt distribution and hence are expected to only 
marginally influence the amount the W boson mass is shifted. 

Our results demonstrate that, for the current level of precision, the approximate cal- 
culation of Ref. |13| is adequate. The small difference in the W boson mass obtained in 
the complete 0(a 3 ) and the approximate calculation, however, cannot be ignored if one 
attempts to measure the W mass with high precision at hadron colliders. This also raises 
the question of how strongly multiple final state photon radiation influences the measured 
W boson mass. So far, only partial calculations exist [p0|1 . A more complete understanding 
of multiple photon radiation is warranted. 

As an alternative to the transverse mass and the lepton pt distribution, the W to Z 
transverse mass ratio, Rm t , has been used recently to extract the mass of the W boson. 
We found that, since the 0(a) corrections to the Z boson transverse mass distribution are 
significantly larger than those to the W Mt distribution, electroweak corrections influence 
Rm t more strongly than the Mt or Pt(£) distribution. 

Finally, we studied how electroweak radiative corrections influence the W cross section, 
the W to Z cross section ratio, Rw/z, and the charge asymmetry of leptons in W decays, 
A(y(£)). As shown in Table [IT], the 0(a) electroweak corrections can reduce the W cross 
section by up to 5% in the presence of cuts. The size of the 0(a) electroweak corrections is 
of the same order as the NNLO QCD corrections [|6|. The 0(a) corrections were found to 
enhance Rw/z by up to 6.5%. QCD corrections, on the other hand, cancel almost perfectly 
in the W to Z cross section ratio and are of (9(1%). The 0(a) electroweak corrections to 
A(y(£)) may not be negligible in view of the projected accuracy of the charge asymmetry in 
future Tevatron runs. 
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TABLE I. Summary of lepton identification requirements. 



electrons 


muons 


combine e and 7 momentum four vectors if 
AR n < 0.2 and if E 1 < 0.15 E e for 0.2 < Ai? e7 < 0.3 


reject events with E~ > 2 GeV 
for AR^ < 0.2 


reject events with E~ > 0.15 E e 
for 0.2 < Ai? e7 < 0.4 


reject events with E y > 6 GeV 
for 0.2 < Ai? M7 < 0.6 



TABLE II. The electroweak if -factor iT w = a°^{pp -> W -> ivX) / 'a Bom (pp -> W -> &/) 
(£ = e, //) and the correction factor to R w/Z , K§ w = RwTz/ R wlz^ with 75 GeV < ™(^ + ^) 
< 105 GeV, for pp collisions at -^/s = 1-8 TeV. Shown are the predictions without and with the 
lepton identification requirements of Table | taken into account. The cuts imposed are listed in 
Eqs. (|36|) and (^). The energy and momentum resolutions used are described in Sec. IIC. 





without lepton id. 


with lepton id. 




requirements 


requirements 


K EW ( pp _^ e + v X) 


0.955 


0.984 


K EW (pp ^+ vX ) 


0.975 


0.947 


Kl W (e) 


1.032 


1.002 


kR w (m) 


1.012 


1.065 
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Born-diagram: 



H {Pi 



(Ji> {Pi 



MPf) 




pure weak contribution: 




virtual 7 contribution: 




real 7 contribution: 

7 




FIG. 1. The Feynman diagrams contributing to W boson production at 0(a 3 ) ($ + : Higgs - 
ghost field, u + ,u y : Faddeev- Popov-ghost fields; the non-photonic contribution to the W self en- 
ergy insertion is symbolized by the shaded loop). An explicit representation of the non-photonic 



contribution to the W self energy insertion can be found in Ref. [26 
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initial: 



final: 



7, k = (1 - z)pi 



h,Ph = zpi 



Si = \(Pi ~ kf - m}\ 
= | - 2k°p° i (l - cos 6) | 
coll. region: < 2(1 — z)(p°) 2 5 fl 



coll. region: 



= + &) 2 — my | 

= |2Jfc°pJ(l - cos 0)| 
< 2(1 - z)z(plf5 e 



FIG. 2. The splitting of the initial and final state fermions i(pi) — > ^(p/i) + 7 and 
^(Pft) ^ f(Pf) + 7 in the collinear region. The hard momentum p^ represents the amount of 
the parent momentum pi j after (before) the emission of a collinear photon. 



p(P) 




FIG. 3. The QED one-loop corrections in deep inelastic lepton - nucleon scattering. 



28 



0.03 



0.02 



0.01 



3 o.oo 



-0.01 



-0.02 



-0.03 



"i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

a ) PP^l +y (7) ISR only 
Vs = 1.8 TeV 



ct(2^3) 

" — 1 



J L 



ct(2->2) + ct(2^3) - CT(Born) 



£j(2-»2) -cj(Born) 



6 a = 0.001 



_l I I I I I I I I I I I I I I I I I L_ 



0.03 



0.02 



0.01 



0.00 



-0.01 



-0.02 



-0.03 



1 1 1 1 1 1 1 1 

- b) pp 


| i i i i | i i i i | i i i i | i i i i 
->l + i>(7) ISR only - 


~_ Vs 


= 1.8 TeV _1 




ct(2^3) 




~ ~1 _ 


L_ ct(2->2) 


+ CT(2^3)-CT(Born) _J 




a(2^2)-a(Born) 


i i i i 1 i i i 


6 S = 0.01 ~ 

1 i i i i 1 i i i i 1 i i i i 1 i i i i 



-3.0 -2.5 -2.0 -1.5 -1.0 -3.50-3.25-3.00-2.75-2.50-2.25-2.00 
l°gio(0 logio(<5 e ) 



FIG. 4. The QED-like initial state corrections to the pp — > l + v(~{), (I = e, fj) cross section 
for yfs = 1.8 TeV as a function of a) S s for 5q = 0.001, and b) 6g for S s = 0.01. Shown are 
a(2 — > 2) — cr(Born), cr(2 — ► 3), and a (2 — > 2) + a (2 — > 3) — <r(Born). The cuts imposed are listed 
in Eqs. ( |36|) and (|37|). The energy and momentum resolutions used are described in the text. 
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FIG. 5. The QED-like final state corrections to the cross section a) a(pp — ► e + iy( , y)) and b) 
a(pp — ► n + u( r y)) for y/s = 1.8 TeV as a function of S s . Shown are cr(2 — > 2) — cr(Born), cr(2 — > 3), 
and cr(2 — ► 2) + a(2 — > 3) — <r(Born). The cuts imposed are listed in Eqs. (|36|) and (|37|). The energy 
and momentum resolutions used are described in the text. 
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FIG. 6. The pp — > e + v('y) cross section for yfs = 1.8 TeV as a function of 5g for 5 S = 0.01 
when electron identification requirements are taken into account and the mass singular terms are 
canceled analytically. Only QED-like final state radiative corrections are included. Shown are 
a(2 -»■ 2) - cr(Born), <r(2 -» 3), and <r(2 -> 2) + a{2 3) - ff(Born). The cuts and lepton 
identification requirements imposed are listed in Eqs. (^) and ([37]), and in Table |. The energy 
and momentum resolutions used are described in the text. 
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50 60 70 80 90 100 25 30 35 40 45 50 



M T (GeV) p T (l) (GeV) 



FIG. 7. Differential cross sections for pp — ► l^v^) at yfs = 1.8 TeV. Shown in part a) is 
the transverse mass distribution. Part b) displays the lepton transverse momentum spectrum. 
The solid (dotted) lines show the distributions for electron (muon) final states including 0(a) 
electroweak corrections. The dashed (dot-dashed) line gives the e + v (/u + z^) Born cross section. 
The cuts imposed are listed in Eqs. (|36|) and (|37|). The energy and momentum resolutions used are 
described in Sec. IIC. The lepton identification requirements of Table | are not taken into account 
here. 
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FIG. 8. Ratio of the C(a 3 ) and lowest order cross sections as a function of the transverse mass 
for a) pp — ► e + u( r y) and b) pp — > ^ + z/(7) at = 1.8 TeV for various individual contributions. The 
upper (lower) solid lines show the result for the QED-like initial (final) state corrections. The upper 
(lower) dotted lines give the cross section ratios if both the QED-like and modified weak initial 
(final) state corrections are included. The dashed lines display the result if only the initial - final 
state interference contributions are included. The cuts imposed are listed in Eqs. ( |36| ) and (|37|). 
The energy and momentum resolutions used are described in Sec. IIC. The lepton identification 
requirements of Table | are not taken into account here. 
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FIG. 9. Ratio of the full 0(a 3 ) and lowest order differential cross sections as a function of the 
transverse mass for a) pp — ► e + u( r y) and b) pp — ► /x + z^(7) at yfs = 1.8 TeV. The dashed (solid) lines 
show the result with (without) the lepton identification requirements of Table | taken into account. 
The cuts imposed are listed in Eqs. @) and @. The energy and momentum resolutions used 
are described in Sec. IIC. 
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FIG. 10. Ratio of the My distributions obtained with the QED-like final state correction part of 
our calculation to the one obtained using the approximation of Refs. [ 13 1 and ]14| for pp — > ^ + z^(7) 
at = 1.8 TeV a) without and b) with lepton identification requirements (see Table |) taken into 
account. The solid and dashed lines give the results for electron and muon final states, respectively. 
The approximate NLO transverse mass distribution does depend on 5 S (see Eq. (|28|)) which is taken 
to be 5 S = 0.01. The cuts imposed are listed in Eqs. (|36[) and (|37|) . The energy and momentum 
resolutions used are described in Sec. IIC. 
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FIG. 11. Ratio of the C(a 3 ) and lowest order W + to Z transverse mass ratio as a function of 
the scaled transverse mass, Xm t , at \fs = 1.8 TeV. The solid (dashed) lines show the result for 
the electron (muon) final state. The ratio without and with lepton identification requirements (see 
Table |) taken into account is shown in part a) and part b) of the figure, respectively. The cuts 
imposed are listed in Eqs. (j36[) and (|37|). The energy and momentum resolutions used are described 
in Sec. IIC. For pp —* £ + l~("f), we in addition require the di-lepton invariant mass to satisfy the 
constraint 75 GeV < m(£ + £~) < 105 GeV. 
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FIG. 12. The charge asymmetry for leptons, A(y(£)), in W — > iv decays for pp collisions 
at y/s = 1.8 TeV. The dashed (dotted) lines show the asymmetry for electron (muon) final states 
including 0(ot) electroweak corrections. The solid line gives the Born prediction of A(y(£)). Except 
for the pseudo-rapidity cut on the charged lepton, the cuts listed in Eqs. (^) and ( |37| ) are imposed. 
The lepton identification requirements of Table | are not taken into account. The energy and 
momentum resolutions used are described in Sec. IIC. 
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FIG. 13. The difference of the 0(a 3 ) and the Born charge asymmetry for electrons (solid) and 
muons (dashed) a) without and b) with the lepton identification requirements of Table | taken into 
account. Except for the pseudo-rapidity cut on the charged lepton, the cuts listed in Eqs. (36) 
and (|37| ) are imposed. The energy and momentum resolutions used are described in Sec. IIC. 
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